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ABSTRACT 

We study the extent to which one can distinguish primordial non-Gaussianity (NG) arising 
from adiabatic and isocurvature perturbations. We make a joint analysis of different NG 
models based on various inflationary scenarios: local-type and equilateral-type NG from adi- 
abatic perturbations and local-type and quadratic-type NG from isocurvature perturbations 
together with a foreground contamination by point sources. We separate the Fisher infor- 
mation of the bispectrum of CM B temperature and polari zation maps by I for the skew 
spectrum estimator introduced bv lMunshi fc Heavens! (|2009f ) to study the scale dependence 
of the signal-to-noise ratio of different NG components and their correlations. We find that 
the adiabatic and the isocurvature modes are strongly correlated, though the phase dif- 
ference of acoustic oscillations helps to distinguish them. The correlation between local- 
and equilateral-type is weak, but the two isocurvature modes are too strongly correlated 
to be discriminated. Point source contamination, to the extent to which it can be regarded 
as white noise, can be almost completely separated from the primordial components for 
/ > 100. Including correlations among the different components, we find that the errors of 



the NG parameters increase by 20-30% for the WMAP 5-year observation, but 
Planck observations. 
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INTRODUCTION 

The statistical properties of fluctuations in the early Uni- 
verse can be used to probe the very earliest stages of its 
history, and provide valuable information on the mechanisms 
which ultimately gave rise to the existence of structure within 
it. This may include evidence for the cosmic inflationary 
expansion. With the recent claim o f a detection of non- 
Gaussianity (|Yadav fc Wandeltl |20Qgh in the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) sky maps, interest in pri- 
mordial non-Gaussianity has obtained a tremendous boost. 

Non-Gaussianity from the simplest inflationary mod- 
els based on a s ingle slowly-rolling sca l ar field is typi- 
cally very small (ISalopek fc Bond! 
Gangui et al.l 1 1994: Acquaviva et al.l 



which the non-Gaussianity of Bardeen's curvature perturba- 

i ; 1|— — I I r 1 

tions is locally characterized (Gangui et al. 1994; Verde et al. 

| 2000l:IWang fc Kamionkowskill^OOOl : iKomatsu fc Spergelll200ll : 

Babich fc Zaldarriagall2004h : 



*(*) = tj>{x) + /nl(^ 2 (x) - (<j> 2 (x))), 



(1) 



1 Il990l: 
2003; 



iFalk et all I l993| : 
Maldacenal hood; 



Bartolo. Matarrese fc Riottoll200a ). Variants of the simple in- 
flationary models can lead to mu ch higher levels of non- 
Gaussianity, such as mult i ple fie lds l|Linde fc Mukhanovlll997l ; 
iLvth. Ungarelli fc W ands 2 Q0j); modulated reheating sce- 
narios jD vali, Gru zinov fc Z aldarriag al 120040 ; warm inflation 
Gupta et alj 120021; iMoss fc Xiond 120071 ): ekpyrotic model 



where cj> is the linear Gaussian part of 4>. This form is mo- 
tivated by the single-field inflation scenarios a nd then many 
mode ls predict non-Gaussianity in terms of /nl (B artolo et alJ 
20041 ) . Optimized estimators of the bispectrum, which is the 
lea ding correlation term in the local form, are introduced 
by iHeavensl l)l998l ) and have been successively developed to 
the point where an estimator for /nl saturates the Cramer- 
Rao b ound for partial sky coverage an d inhomogeneous 
noise l)Komatsu. Spergel fc Wandclt 200E|_ Creminelli et al] 

20061: I Creminelli. Senatore. fc Zaldarriaga] 120071 

Medeiros fc ContaldillgOOd: ICabella et al. 1 12009; biguori et al 



20071 : iKomatsu et al.l 120091 : ISmith. Senatore fc Zaldarriaga 



Kovama et al.ll2007h . 



Different forms are proposed to describe primordial non- 
Gaussianity. Much interest has focused on local-type /nl by 



The local-type /nl is sensitive to the bispectrum with 
squeezed-configuration triangles (k\ <C ~ fe). Several mod- 
els incl uding the inflation scen a rio with non-canonical kineti c 
terms (|Seerv fc Lidsevl 120051 ; IChen. Easther fc Liml |2007h . 
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Dirac -Born-Infeld models ( Alishahiha. Silverstein fc Tond 
l2004l l. and Ghost inflation (|Arkani-Hamed et al.l |2004 ) p~ 
diet large NG signals in equilateral configuration triangles 
(t\ ~ £2 — l 3 ), which is well described wit h equilateral-type 
/nl <|Babich. Creminelli fc Zaldarriagal 120041 ). 

Non-Gaussianity arising from primordial isocurvature (en- 
tropy) perturbations has been discussed in t he context of 
NG field potentials lEnde fc Mukhanovlll997t; , iPeebled Il999l ; 
iBoubekeur fc Lvth 2006; Suvama & Takahashi 2008 ), the cur- 
vaton scenario jlLvth. Unearelli fc Wanddl200 3; Ba rtolo et al.l 
12004 iBeltranl 120081: iMoroi fc TakahashJ I2009T ), modulated 
rehe ating (|Boubekeur fc Creminelli 2006), baryon asymme- 
try (Kawasa ki. Nakavama fc Takahashil 120 091. and the ax- 
ion (jKawasaki et al.ll2008l l. iHikage et"all (f2009l first put ob- 
servational limits on the isocurvature non-Gaussianity using 
WMAP 5-year data. 

In this paper, we make a joint analysis of the different NG 
models to estimate the extent to which one can decode each NG 
information from CMB temperature (T) and E polarization 
(E) maps obtained by WMAP and Planck. We separate Fisher 
information of the CMB bispectrum by different ranges of I to 
study at which angular scale each NG parameter has large S /N 
and correlations among different NG components weaken. This 
idea is based on a new es timator called skew spectrum, which 
iMunshi fc Heavens! <(2009h has introduced to measure a scale 
dependence of NG p arameters, while the commonly-used single 
skewness parameter (jKomatsu. Spergel fc Wandefj2005h gives 
a single value averaged over all scales. The advantage of the 
new estimator is that it retains information on the source of 
the non-Gaussianity, which the commonly-used one does not. 

For our analysis, we adopt a set of cosmological parame- 
ters at the maximum likelihood values for a po wer-law ACDM 
mode l from the WMAP 5-year data only fit l|Dunklev et al.l 
120091 ): n b = 0.0432; Oxim = 0.206; Q A = 0.7508; H = 
72.4 km s _1 Mpc -1 ; r = 0.089; n = 0.961. The amplitude 
of the primordial power spectrum is set to be 2.41 x 10 -9 
at = 0.002Mpc~ 1 . The spectra of isocurvature perturba- 
tions are assumed to be scale-invariant. The radiation trans- 
fer functions for adiabatic and isocurvature perturbations 
are computed using the p ublicly- available CMBFAST code 
(|Seliak fc Zaldarriagal [1999 1. 

This paper is organized as follows: different NG models 
from primordial adiabatic and isocurvature perturbations are 
introduced in $2] Sj3] presents a Fisher matrix analysis of these 
parameters in which we estimate the corresponding error ex- 
pected from WMAP and Planck observations; §4] devotes to a 
summary. 



2 MODELS OF PRIMORDIAL 
NON-GAUSSIANITY 

We consider various forms to describe primordial non- 
Gaussianity from adiabatic and isocurvature perturbations, 
and then provide explicit expressions for the bispectra. 



2.1 Local- Type Adiabatic component 

The bispectrum in the local-typ e NG form (eq.[T]) is written as 
(e.g., iKomatsu fc SpergelllioOll l 



+P (fc 2 )P (fc 3 ) + P (fc 3 )P (fcl)], (2) 

where we rewrite /nl in the equation |T} as /^' ,L ° C . The CMB 
angular bispectra for T, E, and their cross terms are given by 



iAdi,Loc rjrAdi.: 

"xrZ,!i! 2 I 3 — 'JNL 



r 2 dr[/3^ 1 1 (r)/3yf 2 1 (r)a^ 1 (r) 



+Pxt {r)ctY% (r)Pz?s{r) + a*?} (r)/3yf 2 1 (r)/3^, d 3 1 (r)] 



(3) 



where X, Y, and Z denote T or E, and an( A Pxf 1 are 

defined with the adiabatic radiation transfer function g^ 1 as 



OLxi\r) = ~ I k l! dkgxT(k)ji(kr) y 

TV ' 



Pxtir) 



k 2 dkP 4> {k)gxf{k)ji(kr) 



(4) 
(5) 



2.2 Equilateral- Type Adiabatic Component 



The bispectrum in the equilateral-type NG form 
is characterized by the NG parameter /**' Eq 
l|Babich. Creminelli fc Za ldarriaga 2004) as follows: 

S Adi ' Eq (fci,fe2,fe 3 ) = 6/ N ' di ' Eq [-P (fc 1 )P (fc 2 ) - P^(fc 2 )P^(fc 3 ) 
-P (fc 3 )P (fci) - 2{P (fci)P (fc 2 )P (fc 3 )} 2/3 



+{[p*(fci)p*(fa) 3 p*(fcs 



,3-11/3 



+ (5 perm.)}]. 



The CMB angular bispectra in this form are given by 



, Adi,Eq 
u XYZ,l 1 l 2 l 3 



6 /nl 



Adi,Eq 



2 dr[-P^l(r)P^(r)a^(r) 



/oAdi / \ Adi / \^>Adi/ \ Adi/ \ /.Adi/ \/->Adi/ \ 

-Pxh {r)a Y i 2 {r)Pzi 3 \r) - a xh (r)(3 Y i 2 (r)Pzi 3 (r) 

-20r)#£(r)«££(r) 
+{/3^(r)7rf 2 i (r)4f 3 1 W + (Sperm.)}], (7) 



where 



7 £f(r) = ±J fc 2 dfcP 1/3 (fc) 5 £f(fc)j ; (fcr), 
6 x f(r) = - /fc 2 dfcP 2/3 (%£f(fc)j ; (fcr). 



(8) 
(9) 



2.3 Isocurvature Components 

Here we consider an isocurvature perturbation S between 
axion-type cold dark matter (CDM) and radiation, which is 
uncorrelated with adiabatic perturbations, defined as 

SpcDM 3<5p 7 



S = 



PCDM 



4 Pl 



(10) 



where pcdm is the CDM energy density and p 1 is the radiation 
energy density. The fractional isocurvature perturbation fs is 
defined as 



fa 



Ps(fco) 



P((ko) + Ps(ko) 



(11) 



B- 



c (k 1 ,k 2 ,k 3 )=2f^° c [P4k 1 )P c ,(k 2 ) 



where P^ and Ps represent the power spectra of £ and S and 
fco is set to be 0.002Mpc _1 . At linear order, $ (eq. [T]) is re- 
lated to £ by <3? = (3/5)£. The definition of fs is same as the 
commonly used parameter a |Bean et al.l|2006l ). The current 
observational limit on fs is 0.067 (95% CL) for t he axion-type 
isocurvature perturbation (|Komatsu et alj|2009h . 
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2.3.1 Local-Type Isocurvature Component 



2.4 Point Source Component 



We consider two different forms for isocurvature non- 
Gaussianity. One is the same local form as the adiabatic one 
(eq.0): 



Unmasked point sources (e.g., radio galaxies) generates an ad- 
ditional non-Gaussianity in observed CMB maps. Assuming 
them to be Poisson distribution, b x s YZ ^ ll2l3 ^ s a constant. 



B^'^fci.fe.fca) = 2/^ Loc [P,(fc 1 )P I) (fc 2 ) 

+P v {k 2 )Pr,{k 3 ) + P v (k 3 )P v (k 1 )], (12) 

where r\ denotes the Gaussian part of S with the amplitude 
of Pr/(k) is normalized bv fa (ec\.\TT\). The pa rameter /J," ,Loc 
corresponds to Q 2 /n|° in iHikage et alj (|2009t ). We obtain the 
CMB bispectrum as 



h XYZ,i^H = 2 /nl / r dr[f3 xh (r)p Yl2 (r)a zh (r) 

+Px° h {r)a% (r)/3gf 3 (r) + a& M/3?f 2 (r)/£& (r)] , (13) 
where a^? and /3^° 

are defined with the isocurvature radiation 
transfer function gx° as 



7T ' 



dkP v (k)gx°(k)ji(kr) 



(14) 
(15) 



3 FISHER INFORMATION ANALYSIS FOR 
SKEW SPECTRUM 

We make Fisher information analysis of the different NG com- 
ponents introduced in the previous section to estimate the error 
expected from WMAP, Planck and noiseless ideal observations. 

The Fisher matrix for the CMB bis pectrum in the weakly 
non-Gaussian, all-sky limit is written as l|Babich fc Zaldarriagal 
12004 lYadav. Komatsu.fe Wandeldl2007M 



F K 



l 

E 

2<h<h<l 



(21) 



XYZ PQR 



Xb 



(f~t -1\XY Z\PQR,j (r,r,\ 



where i and j denote each NG component and the factor -Z; 1 ; 2 ; 3 
is defined as 



2.3.2 Quadratic-Type Isocurvature Component 

When the linear Gaussian term is negligible compared with the 
quad ratic term, the isocurvatu re perturbation has a x 2 from 
(e.g., iLinde fc Mukhanov|[l997h : 



S = a 2 - (a 2 ), 



(16) 



where a obeys Gaussian st atistics. This form has been stud- 
ied in the context of axion ( Kawasaki et al. Il2008l ) and curva- 
ton scenarios (ILanglois. Vernizzi fc Wandsl 12008). The bispec- 
tra are calculated as l|Komatsull2002l ) 



(P) 



r>Iso,Quad/, i i \ 8 f rf 3 p 

B (k 1 ,k 2 ,k 3 ) = -J L _ i ^P a 

x [P CT (|ki + p|)P CT (|k 2 - p|) + P CT (|k 2 + P |)P CT (|k 3 - p|) 

+P CT ([k 3 +p|)P CT (]ki-p|)], (17) 

where a finite box-size Lbox gives an infrared cutoff. To avoid 
assumptions at scales far beyond the present horizon Hq 1 , we 
set I/box = 30Gpc. Th e equation (|17|) is approximately given 
bv lHikage et all (|2009f) as 

^7/1*3= 2 / r 2 dr[/3^^(r)$Ur)aX(r) 

+Px° h (r)ay-l (0/3zT' Quad (0 + <*% (r)ffi' Qnad (r)$? 3 (r)], 



(18) 



where 



/3j°' Quad (r) 



Px°(r) 



k a dkP a {k)g%l(k)ji(kr) 



(19) 
(20) 



Non-Gaussianity in this form is characterized by fs feo.|ll|) 



(2Zi + l)(2Z 2 + l)(2Z 3 + l) / h h h 



4tt 







(23) 



The sums over XYZ and PQR are just TTT when using CMB 
temperature maps only (T only), but are eight combinations 
{TTT, TTE, TET, ETT, TEE, ETE, EET, EEE) when both 
CMB temperature and E polarization maps are used (T&E). 
The Fisher matrix at each Z, F ; !J , is associated with the skew 
spectrum estimator for th e i-th NG component, Si, defined as 
|Munshi fc Heavensl l2009) 



SI = 



21 - 



1 ^ 

2<l 1 <l 2 <l 



XYZ PQR 



,i tr ~, -1^XYZ\PQR U obs 

Xt>XYZ,lil2l\^ OV )hl 2 l OpQR.hhh 



(24) 



where b 



denotes the observed bispectrum. The 
relation to the single skew ness estimator 5p rim 
<|Komatsu. Spergel fc Wandeldl2005h is 



5>j+i)s? 



(25) 



When non-Gaussianity is small, the covariance matrix is ap- 
proximately given by 



l~, XY Z\PQR a n XP n YQ n ZR 



(26) 



where A !l!2i3 is 6 (Zi = h = h), 2 (h = h, h = h, or h = h), 
and 1 (Zi 7^ h 7^ Z3) and C^ Y represents the CMB power 
spectrum from purely adiabatic perturbations including ob- 
servational noise AT; : 



Cf Y = - / k 2 dkP^k) gxl \k)g$t(k)+Nl 



(27) 



We consider three different noise/beam functions: an ideal 
case without noise/beam ("Ideal"); WMAP 5-year V+W band 
coadded map ("WMAP5"); Planck's expectations after two 
full sky surveys for 14 months ( "Planck" ) using all of nine fre- 
quency channels. Noise is assumed to be homogeneous white 
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Table 1. Diagonal components of the Fisher matrix F" summed 
up to I = 2500. The different noise/beam for WMAP5, Planck, and 
Ideal are considered. 





Adi,Loc 


Adi,Eq 


Iso,Loc 


Iso,Quad 


WMAP5 (T only) 


2.7x10" 


-3 


7.2x10- 


-S 


2.9x10" 


-4 


6.8x10 s 


Planck (T only) 


3.7x10" 


-2 


2.3x10- 


-4 


3.1x10" 


-4 


7.6xl0 3 


Ideal (T only) 


9.0x10" 


-2 


3.4x10- 


-4 


3.1x10" 


-4 


7.7xl0 3 


WMAP5 (T&E) 


3.0x10" 


3 


7.5x10" 


-S 


3.1x10" 


-4 


7.2xl0 3 


Planck (T&E) 


5.8x10" 


2 


4.8x10- 


-4 


8.9x10" 


-4 


2.6xl0 4 


Ideal (T&E) 


3.6x10" 


1 


2.9x10" 


-3 


3.9x10" 


-3 


1.3xl0 5 



noise and N* Y — when Noise/beam is coadded at 

each I with the inverse weight of the noise variance in each 
frequency band or differential assembly. Planck's noise/beam 
information is obtained from http://www.rssd.esa.int/Planck). 
The fraction of sky / s k y is set to be 1 in this analysis. 

Figure [1] shows the diagonal component of the Fisher ma- 
trix Ff 1 (eq. [22]). It represents the square of signal-to-noise 
ratio (S/N) 2 for i-th NG component at I without correlations 
among different NG components. The adiabatic components 
have an increasing trend of S/N at higher I. The majority of 
the signal of the isocurvature components in temperature maps 
come from the large-angular scale (Z < 100), where isocurva- 
ture perturbations produce larger CMB fluctuations than adia- 
batic perturbations. A phase difference in acoustic oscillations 
between adiabatic and isocurvature modes provides a distinct 
signature seen around Z ~ 300, which is important particularly 
when polarization maps are included. Table [1] lists the values 
of the diagonal components of the Fisher matrix summed over 
Z up to 2500, at which Planck estimates are enough saturated. 

Figure [5] shows the cross-correlation coefficient r; defined 
as Ti = I '(F" Fi') 1 ' 2 . The local-type adiabatic and isocur- 
vature components are strongly correlated, but the phase dif- 
ference of acoustic oscillations weakens the correlation, as seen 
especially around I ~ 200. The correlation between the local- 
type and the equilateral-type components becomes weak at 
I > 100. The two isocurvature components are almost com- 
pletely correlated over all scales. The correlation with the point 
source component is very weak for I > 100. 

Figure [3] shows la error contours (Cramer-Rao bound) 
for a pair of NG parameters for WMAP5 (T only), Planck (T 
only), and Planck (T&E). The errors expected from WMAP5 
(T&E) is almost same as those from WMAP5 (T only). The 
rest of NG parameters other than two plotted are fixed to be 
zero. The local-type adiabatic and isocurvature components 
are correlated with the correlation coefficient r = 0.43 for 
WMAP5, r = 0.23 for Planck T only, r = 0.20 for Planck 
T&E when Z is summed up to 2500. We see that the local- 
type and the quadratic-type scale-invariant isocurvature com- 
ponents are difficult to be differentiated even using Planck 
data. The local-type and the equilateral-type adiabatic com- 
ponents are weakly correlated (r = 0.12 for WMAP5, r = 0.17 
for Planck T only, r = 0.22 for Planck TfcE), which is consis- 
tent w ith the previous work (Babi clu Creminelli fc Zaldarriagal 
2001). The point source component is almost uncorrelated with 
the other primordial components (r < 0.08 for WMAP5 and 
r < 0.03 for Planck), which is consistent with the previous 
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Figure 1. Diagonal components of the Fisher matrix F" (eq. [22]). 
From top to bottom, the local-type adiabatic (Adi,Loc), the 
equilateral-type adiabatic (Adi,Eq), the local-type isocurvature 
(Iso,Loc), and the quadratic-type isocurvature (Iso,Quad) compo- 
nents are plotted. Left panels are for T map only, but right panels 
are for T&E maps. Noise/Beam is for Ideal (solid), WMAP5 (dot- 
ted) and Planck observations (dashed). 



work ([Komatsu fc Spereell l200lh . Table [2] lists the errors of 
the NG parameters without and with correlations among all 
of other parameters except for the quadratic-type isocurvature 
component. Polarization maps are found to be very important 
to constrain the isocurvature NG as well as adiabatic NG. The 
increase of the errors due to the correlations mainly between 
adiabatic and isocurvature modes is 20-30% for WMAP5, but 
less than 5% for Planck observations. 
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Figure 2. Cross-correlation coefficients rV = F^ 3 /(F"F 33 J 1 / 2 where i and j denote the local-type adiabatic (Adi,Loc), the local-type isocur- 
vature (Iso,Loc), the equilateral-type adiabatic (Adi,Eq), the quadratic-type isocurvature (Iso,Quad), and the point source (PS) components. 
Left of each panel is temperature T only; right is T+E, including E polarization. Noise/Beam is for Ideal (solid), WMAP5 (dotted) and Planck 
observations (dashed). 
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Figure 3. la error contours of a pair of NG parameters expected from WMAP5 T only (solid circles), Planck T only (dotted circles), and 
Planck T&E (dashed circles). The rest of NG parameters other than two plotted are fixed to be zero. 



Table 2. la errors of each NG parameter from the joint analysis 
of all of the NG components except for the quadratic-type isocurva- 
ture components. The values without parentheses denote the limits 
with other parameters fixed, but those with parentheses denote the 
limits including the correlations such that the other parameters are 
marginalized. 





WMAP5 


Planck 




T only 


T only 


T&E 


a rAdi,Loc 

fl % r 

a /-lso,Loc 

* NL 
A # Adi,Eq 

J NL 

A6 p S[10" 28 ] 


19 (23) 
59 (82) 
117 (149) 
438 (441) 


5.2 (5.5) 
57 (62) 
66 (71) 
0.13 (0.13) 


4.1 (4.3) 
34 (35) 
46 (48) 
0.024 (0.024) 



4 SUMMARY 

We have presented a detailed analysis of the possibility of ex- 
tracting information about non-Gaussianity from various infla- 
tionary models. We consider four different-type primordial NG 



models: local-type adiabatic, equilateral-type adiabatic, local- 
type isocurvature, and quadratic-type isocurvature models to- 
gether with point source contamination. The adiabatic and 
the isocurvature modes are correlated, but the difference in 
the phase of the corresponding acoustic oscillations breaks the 
degeneracy. The local-type and quadratic-type scale-invariant 
isocurvature components are difficult to separate even using 
Planck data. The correlation between the local-type and the 
equilateral-type adiabatic modes is weak. The point source 
(white noise) contamination does not pose a threat as it is 
uncorrelated with any of the /nl parameters, although a high- 
resolution experiment will be more suited to get rid of such 
contamination. Our results are based on noise models from 
WMAP and Planck and we compare them to ideal noise-free 
and all-sky reference observations. The increase of the error for 
the non-Gaussian parameters due to the correlations is 20-30% 
for WMAP5 and 5% for Planck. 

Secondary anisotropics other than point sources can 
contaminate the estimation of primordial non-Gaussianity. 
The cross-contamination of various inflationary contributions 
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against secondaries such as Sunyaev-Zeldovich effect (SZ) or 
Integrated Sachs- Wolfe effect (ISW) which are potentially ob- 
servable with Planck data will be present elsewhere. 
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